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Two-particle wave function in four dimensional Ising model 

T. Yamazaki a 

a RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973, U.S.A. 

An exploratory study of two-particle wave function is carried out with a four dimensional simple model. The 
wave functions not only for two-particle ground and first excited states but also for an unstable state are calculated 
from three- and four-point functions using the diagonalization method suggested by Liischer and Wolff. The 
scattering phase shift is evaluated from these wave functions. 
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1. Introduction 

Calculation of the scattering phase shift is im- 
portant to understand scatterings and decays of 
hadrons. The phase shift in I = 2 ttt: scattering 
system [3] was evaluated with the finite volume 
method |1I2| derived by Liischer. In the finite 
volume method two-particle wave function plays 
an important role, because the phase shift is ex- 
tracted by the wave function. 

Study of the wave function was carried out by 
Balog et al. @|. Using two-dimensional statisti- 
cal model they evaluated the phase shift from the 
wave function for several two-particle states. Re- 
cently Ishizuka et al. and CP-PACS collaboration 
studied the wave function in I = 2 -kit scattering 
system They extracted the scattering length 
and scattering effective potential from the wave 
function for the ground state. Here calculation of 
the wave function for the first excited and an un- 
stable states is attempted with a four dimensional 
simple Ising model. It is also aimed to evaluate 
the phase shift from these wave functions. 

2. Wave function 

Liischer proved that the wave function <j>(r) sat- 
isfies effective Shrodinger equation, 

(V 2 +p 2 )0(r) = J d 3 rU E (r,r')(f>(r'), (1) 

where r is relative coordinate of two particles, p 
is the relative momentum, and Ue (r, r ) is the 
Fourier transform of the modified Bethe-Salpeter 
kernel introduced in ref. pQ. The effective poten- 
tial UE{r,r ) depends on energy E and decays 
exponentially in r and r . The wave function sat- 



isfies the Helmholtz equation (V 2 + p 2 )4>(r) = 
in r > R. The R is the effective range where 
Ue(t,T ) becomes sufficiently small in exterior re- 
gion of R. Liischer found general solution of the 
Helmholtz equation in a finite volume, 
G(r) = (l/L^e^/(k 2 -p 2 ), (2) 

k 

where k = (2w/L) ■ n and n is an integer vector. 
We can extract p 2 by fitting 4>(r) in r > R with 
G(r), and can then evaluate the phase shift with 
the finite volume method. 

3. Methods 

A simple model, which is constructed with a 
lighter mass field 7r coupled to a heavier mass field 
a with three-point coupling, is employed to cal- 
culate an unstable state. This model has been 
successfully used to observe a resonance by Gat- 
tringcr and Lang |B] in two dimensions, and by 
Rummukainen and Gottlieb [7j in four dimen- 
sions. 

Three- and four-point functions, Fi(f,t) = 
(0\W(f)[Oi(t) - O t {t + 1)]|0) and C ZJ (t) = 
(0\O t (0)[Oj(t) - Oj(t + 1)]|0) for i,j = 0,1 and 
a, are calculated to obtain the wave function for 
some states. The operator Oi(t) is 7T7t operator 
with p — 0, (2ir/L) (i = 0,1) and a operator 
(i = a), and the subtraction [Oi(t) — Oi(t + 1)] is 
performed to eliminate the vacuum contribution. 
The wave function operator W(r) is defined by 
W(r) = (l/L 3 )^^ R 7r(R[rl + A ? )vr(A), and the 
R is an element of cubic group, and summation 
over X and R projects to A + sector. The lattice 
size is L 3 x T = 20 3 x 64 and the number of the 
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Figure 1. State overlaps of operators at several 
t. 

configurations is 405 xfO 3 . In this work the en- 
ergy of the first excited state is larger than that 
of the a state. 

The wave function is defined by (j) a (r) — 
(0\W(r)\a) for the a = 0,1 and a state. Using 
4> a (r) and the state overlap of the operator V a i — 
(a\Oi\0), one has Ft(?,t) = Ea^(0 A «W^, 
where A a (t) = (1 - e- B °)e~ E °' + (1 - 
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-E a {T-t)_ is also easy to show that 

Cij{t) = Ea^r^«(*)^»j- Tne assumption of 
this calculation is that higher energy states than 
the first excited state do not contribute to Cy (t) 
and Fj(f, t). To obtain the wave function, at 
first the overlaps are extracted by the diagonal- 
ization method jS] with the correlation function 
matrix C- 1/2 (t )C(t)C~ 1/2 (t Q ), where t is ref- 
erence point. The diagonalization also yields the 
energy E a . Then the wave function is extracted 
by the projection 

^=[^mt)Vr 1 ]/\ZF i (f ,t)V^ 1 \ (3) 

i i 

apart from the normalization, where the normal- 
ization point ro is chosen ro = 1. 

4. Results 

Fig-Qdisplays the overlaps determined at some 
t. The i — operator is almost dominated by the 
ground state, while other operators has contribu- 
tions from each state. All the overlaps are stable 
at small t region, so that Vi a is determined at 
t = 1. Using the overlaps it is possible to cal- 
culate the wave function not only for the first 
excited state but also for the a state as shown 




Figure 2. Wave function calculated from eq. (3) 
for ground, first excited and a states. i*i(r) cor- 
responding to each state is also presented. 

in Fig. The projection eq. (3) is carried out 
with Fi(r,t) at t = 3 for a = 0, 1, and at t = 4 
for a = a. The function Fi(r,t) normalized at ro 
is also displayed in the figure for comparing ones 
before and after the projection. The difference 
between <p a (r) and Fi (r, t) for the ground state 
is small, while those for the first excited and a 
states are large, as expected from the operator 
overlaps in Fig. ^ 

The effective range is required for fitting of 
4> a (r) in r > R using eq. (2). In Ref. [S] the 
effective range is estimated from the quantity 
CV 2 (l) a {r))/(p a {r), which is approximately effec- 
tive potential. However, the effective range can- 
not be estimated from the quantity, because the 
statistical noise of the quantity is very large. 
Hence the fit of (f) a (f) is carried out by assuming 
the effective range R = 7,9 and 8 for a = 0, 1 
and a states. The effective range depends on 
the state, because Ue{t, r ) in eq. (1) depends on 
the energy. The fitting parameters are an overall 
constant of eq. (2) and the relative momentum 
p 2 . Fig. 01 shows the fitting result G a (r) for each 
state. The fitting results for all states are con- 
sistent with the wave functions. This figure illus- 
trates that the wave function in r > R for the 




Figure 3. Fitting result of wave functions for 
ground, first excited and a states. 

a state can be described by the general solution 
of the Helmholtz equation. The consistency of 
the effective ranges is checked by the deviations 
of 4>a.{r) from G a (r). As expected, the deviation 
vanishes in r > R. 

As shown in Fig. the phase shift is evalu- 
ated using p 2 determined from the wave function. 
In the figure the phase shift calculated with p 2 
obtained from the energy E a — 2^f m 2 + p 2 is 
also presented. All the phase shifts are consis- 
tent with the two methods. The results obtained 
from the wave function for the ground and hrst 
excited states have smaller error than those ob- 
tained from the energies. This feature is also seen 
in refs. 0] and The wave function result for 
the a state, however, has larger error. The reason 
is not well understood at present. In order to un- 
derstand this feature, more detailed investigation 
is needed for the wave function of the a state. 

5. Conclusions 

Calculation of the two-particle wave function 
for the first excited and an unstable states is 
demonstrated with the diagonalization method. 
It is found that the wave function in r > R for 
an unstable state can be described by the gen- 
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Figure 4. Scattering phase shift (degrees) ob- 
tained from wave function (W) and energy (E) 
for ground, first excited and a states. 

eral solution of the Helmholtz equation as same 
as that for two-particle states, and the phase shift 
is extracted from these wave functions. 

The numerical calculations have been carried 
out on workstations at Center for Computational 
Sciences, University of Tsukuba. 

REFERENCES 

1. M. Liischer, Commun. Math. Phys. 105 
(1986) 153. 

2. M. Liischer, Nucl. Phys. B354 (1991) 531. 

3. CP-PACS Collaboration: S. Aoki et ai, 
Phys. Rev. D67 (2003) 014502; CP- 
PACS Collaboration: T. Yamazaki et ai, 
Nucl. Phys. BfProc. Suppl.)129 (2004) 191; 
|hep-lat7 0402025. 

4. J. Balog et.al, Phys. Rev. D60 (1999) 094508; 
Nucl. Phys. B618 (2001) 315. 

5. N. Ishizuka and T. Yamazaki, Nucl. Phys. 
B(Proc. Suppl.)129 (2004) 233; CP-PACS 
Collaboration: S. Aoki et ai, these proceed- 
ings. 

6. C. R. Gattringer and C. B. Lang, Phys. Let. 
B274 (1992) 95; Nucl. Phys. B391 (1993) 
463. 

7. K. Rummukainen and S. Gottlieb, Nucl. 
Phys. B450 (1995) 397. 

8. M. Liischer and U. Wolff, Nucl. Phys. B339 
(1990) 222. 



